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Irreversibility and acausality of a sub-system are established in exactly soluble harmonic models 
with reversible and causal dynamics. It is shown that initial conditions, imposed on some dynamical 
degrees of freedom may break time reversal invariance for other degrees of freedom. This happens if 
observations carried out in any large but finite amount of time can not resolve the spectrum of the 
eliminated degrees of freedom, namely when the spectrum has a condensation point at the ground 
state. Acausality follows due to the dominance of the dynamics by almost time-independent modes. 
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I. INTRODUCTION 



o 

All fundamental forces apart of weak interactions are time reversal invariant. Hence we expect that effective theories, 
O ' constructed for energies well below the energy scale of weak interactions preserve the symmetry with respect to time 
reversal. But coarse graining even of a reversible dynamical system [1] already leads to a definite time arrow [2], 
a dynamical breakdown of the time reversal invariance. Irreversibility presents itself usually by the mixing of time 
, derivatives of odd and even powers in the linearized effective equations of motion. These terms induce other changes 
in the dynamics, for instance they may generate complex poles for the Green-functions in energy and endanger the 
r— I \ expected causal structure and the unitarity of time evolution. Such a relation between dissipative forces and acausality 
'"pi ' is well know in case of the Abraham-Lorentz-Dirac radiation reaction force in classical electrodynamics [3] . 
' ] The simplest phenomenological treatment of dissipative forces is the insertion of a term into Newton equation which 
is linear in the velocity. The resulting equation violates time reversal invariance but remains causal. It is difficult even 
(-H ' to imagine how dissipation may lead to acausality since one would think that we can always construct the solution of 
I— ' , the canonical equations of motion by a simple integration in time. A loophole appears in this argument when applied 
to a system involving infinitely many soft degrees of freedom. Weaker restoring force to the stationary position means 
slower motion and we may need infinitely long time to get all degrees of freedom moving. We encounter two limits in 
' this procedure, the long observation time T — oo and the large system N ^ oo which may not commute, rendering 
' the time evolution non-unique and requiring the consideration of certain details of the observation to decide the order. 
Cn| The strategy of effective theories will be used in the present work to trace the origin of acausal and irreversible 

C^l ' effects starting with an underlying reversible, unitary dynamics. The separation of the degrees of freedom into a 
, system and its environment and the elimination of the environment variables by means of their equations of motion 
allows us to define, at least formally, the effective dynamics of the system. It is found that irreversibility and acausality 
may appear when infinitely many eliminated environmental degrees of freedom remain unresolved by the observation. 
These degrees of freedom provide a sink for the dissipated energy and a mechanism to smear the impact of external 
perturbation in time and thereby generate irreversibility and acausality. 

In particular, we consider exactly soluble systems where the Lagrangian is quadratic in the coordinates and the 
velocities but we believe that our qualitative conclusions remain valid for interactive systems, too. The possibility of 
a phase transition opens here by coupling an observable to infinitely many normal modes. The way this infinity is 
reached may generate singular dependence on the parameters of the model, the hallmark of phase transitions. Such a 
broad framework has already been used to explore dissipative dynamics [4, 5]. Up to our best knowledge the resulting 
acausality has not been commented or traced back to its origin. The generic harmonic model consists of a set of 
linearly coupled harmonic oscillators and can be characterized by its spectral functions. A more restricted class of 
models is where the environment is realized by a free field, coupled linearly to a harmonic oscillator, considered as 
"system" . In such a model the spectral function of the environment is determined by the space-time symmetries of the 
field theory in question. We shall use this model together with the generic one given in terms of harmonic oscillators 
and analyze both of them on classical and quantum level. The causal structure is naturally identical for quantum and 
classical harmonic systems but the issue of unitarity can be considered in quantum mechanics only. 

Reversibility and causality are always found to be intact for systems with discrete excitation spectrum, as expected. 
But time is needed for the observations to resolve the discreteness of an excitation spectrum. When the time available 
for observation is not sufficient then the discrete spectrum appears to be spread. This cause a drastic change if there 
is a condensation point of the spectrum at the ground state because no observation carried out in finite amount of 
time can resolve this part of the spectrum. It is the contribution if these unresolved modes which makes the effective 
system dynamics irreversible and acausal, confirming the important role of information loss in dissipative systems. 
The first part of the paper builds up the effective theory for the system by eliminating the environment. The 
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classical models consisting of harmonic oscillators or a scalar field as environment are introduced in Section II. The 
quantum! case is considered in Section III starting with a brief summary of the Closed Time Path (CTP) formalism 
[6] which is needed to preserve the transparency offered by Wick theorem, a necessary tool to treat realistic models. 
Physical considerations are presented in the second part, starting with the interpretation of the irreversibility as 
radiation damping in Section IV. The similarity with spontaneous symmetry breaking is addressed in Section V. 
Finally, the conclusions are summarized in Section VI. 



II. CLASSICAL MODELS 



We introduce in this section two simple models consisting of linearly coupled harmonic oscillators. One of them 
is called system and the rest represents the environment. The spectrum for the latter can freely be chosen in the 
first model and is restricted by space-time symmetries in the second case. A further notational difference is that 
the spectrum is kept discrete in the first case for the sake of simplicity and the second, field theory based model is 
considered in infinite space with continuous spectrum. 



A. Harmonic oscillators 



The simplest nontrivial quadratic model consists of a system of linearly coupled harmonic oscillators [5]. The 
Lagrangian is written as 
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(1) 



where y and x„ denote the system and the n-th harmonic oscillator coordinate, respectively. The source j is introduced 
to diagnose the system dynamics. The equations of motion 



y = j- rmoly - ^ gnXn, 



my 

mXn = -mu;^Xn - QnV, 

can easily be solved by successive elimination. We impose the initial conditions 

y{ti) = y{ti) = 0, 



(2) 



(3) 



for the system and x„(ii) = Xhniti), Xn{ti) = XhniU) with Xhn{t) = SR^nC*'^'** for the environment at the initial time 
t = ti and go into the limit ti = — oo to recover continuous frequency spectrum in the Green-functions. The solution 
of the equations of motion for the Fourier transform 



.) = / 



dxe^*f{t) 



of the environment coordinates is 

a;„(w) = 



The system coordinate is expressed in terms of the retarded Green-function 



as 



roo 

y{t)= / dt'Cit-t' 



^gnXhnit') - j{t') 



where the Fourier transform 



Go^(w + «e)-I]'^(w) 



(4) 



(5) 



(6) 



(7) 



(8) 
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is obtained by means of the free inverse propagator 

G^\u) = m{i^^-wl) (9) 

and the self energy 

E» = l^^— 4 2- (10) 

m ^ (u) + ley — 

Each order of the geometric series resulting from the expansion of the propagator (8), 

G'{t,t') = G'o{t,t')+ J dhdHGl{tMWQ{U,t2)Gl{t2,t') + --- (11) 

is non- vanishing for t > t' only and G^ is causal as long as the expansion converges. This is the case when the G"'(a;) 
has a unique analytic extension in the complex frequency plane, in case of discrete spectrum. Hence causality is assured 
for finite systems or for an infinite set of harmonic oscillators with spectrum without condensation point. When the 
spectrum possesses condensation point then the frequency Fourier integrals of G^{ui) with the self energy (10) in the 
denominator may not have uniform convergence causing a dependence on the order the frequency integration and 
the mode summation are carried out. When the frequency integral is made first and the summation over the; modes 
is performed next then the pole structure of (8)- (10) assures causality which may be lost by following the opposite 
order. 

An alternative way to find the retarded Green-function is based on the normal modes 

*J=E^7>" (12) 

n 

which diagonahze the Lagrangian, 

m^-—L:ilj+jJ2AojXj. (13) 
The retarded Green-function for the system coordinate is of the form 

^2 

G'-{uj) = y— ^ ^ (14) 



with Y,- Ah = 1 yielding 



G^{t) = -Q{t)J2^'^^. (15) 

j ■' 

The causal nature of this Green-function is again due to the fact that the frequency integral was dealt with before 
the summation over j in this equation. 

B. Spectral strengths 
It is advantageous at this point to introduce some spectral functions. The spectral strength 

Pem = E^i^-n)^ (16) 
n 

of the environment allows us to write the self energy as 

-'M = /<'0^r».. (.7) 



(w + ief - 02 • 
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Whenever this self energy is used we tacitly commit ourselves to a normal modes summation carried out before 
frequency integration. 

The spectral function of the normal modes is defined by 

p^{n)=Y^S{il-u:^)^, (18) 



3 ■' 



and the Green-function 



is naturally retarded, 

G''{t) = -e{t) j dfl2pm{^)smnt, (20) 

when the frequency integral is made before the integration over the spectral variable. 

To find a model for dissipativc phenomena it is usually assumed that the environment spectrum has a condensation 
point at w = 0. The long time phenomenology of the resulting system-environment interactions depends on the 
asymptotic of the spectral function at vanishing frequency, assumed to be pd^^) = 0(nP). An environment with 
p = I called Ohmic and it separates the sub- and super-Ohmic regimes with p < 1 and p > 1, respectively [5]. The 
spectral function should approach zero for large frequencies, lima^aa Pe{^) = to avoid UV divergences. A simple 
realization of an Ohmic environment is given by the Deby-suppressed Drude model, 

where ge and rUe are averaged coupling strength and mass parameter of the environment. The self energy (17) is now 



and Green-function (8) 



S'-(a;) = . (22) 



G^{io) = (23) 



turns out acausal in general. 

One may construct similar phenomenological spectral function for the normal modes, as well. We assume that the 
condensation point of the normal mode spectrum is located at the lowest frequency u^mi^ '^o) and write 

p^iSi) = &in - ^m) ~.~^f.^ (24) 

mi2[S2^ -I- (12 — cJm) \ 

The Green-function can easily be calculated for Um = and we find poles for both sides of the real frequency axis in 

The generic, formal feature of a dissipative model is clearly shown by these examples. We start with a "microscopic" 
model with a large but finite number of degrees of freedom which is a closed system with time reversal invariant and 
causal dynamics. These latter properties might be lost for an infinite system. It is usually the thermodynamical 
limit which produces continuous spectrum and the emerging the condensation points in the spectrum may generate 
complex poles for the Green-function leading to a breakdown of time reversal invariance. When some of these poles 
end up on the "wrong" side of the complex frequency plane then acausality is observed. 



C. Classical field as environment 



To place the model into a more realistic context we consider spectral strength which is fixed by symmetries. The 
harmonic oscillators of the environment will then be described by a free real scalar field (/>(f, x) in space dimension d 
and of mass /x. The system coordinate y is coupled to the field at the origin and the action is chosen to be 



S[cp,y]= / dt 



y''{t)-gy{t)cj>{t,0)+j{t)y{t) 



dtd'^2 



■ (25) 
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The coupling constant will be written as 17 = yfmujQl'^l'^ ^ in terms of the parameter I of dimension of length. The 
equation of motion for the environment normal modes, 

^{t,k) = J d^xe-'''''(l){t,x) (26) 

is 

[d^+io^{kMt,k) = -gc'y{t) (27) 

where w(fc) = c^J ^^0- + fe^. The corresponding initial conditions are chosen to be (/>(fj, fe) = (phiU, k) and dt4){ti, k) = 
dtft^hiti, k). The environment normal coordinates are inserted into the system equation of motion 

^^{t,k)+j{t) (28) 

yielding the solution 

/•oo 

y{t)= / dt'G''{t-t')[gMt',x = 0)-jit')] (29) 



which satisfies the initial conditions y{ti) = yifi) = 0. The Fourier transform of the is given by Eq. (8) except 
that the self energy is now given by 

^>^=y (2^(j7w^^- (^0) 

We introduce the spectral weight 

which allows us to write the self energy as in Eq. (17). In the case of gapless environment, = 0, we have 

p,{n) = ^(-Wa (32) 



27r V c, 

with pi = P2 = \i p3 = l/27r, and = l/(27r)^. The thcrmodynamical limit makes the spectrum continues and the 
number of soft modes is just right to make up an Ohmic dissipation in 3 spatial dimensions. 



III. QUANTUM MODELS 



After having seen the mathematical possibility of generating irreversible or acausal dynamics wc turn to quantum 
systems. The retarded solution was manufactured in the classical case and the eventual acausality in the retarded 
Green-function is to be considered by suspicion. On the contrary, the solution of the Schrodinger equation will 
be obtained by functional method, applicable for either finite or infinite system and the retarded Green-function is 
thereby constructed in a systematic, unique manner. 



A. CTP scheme 



The widely employed formalism of quantum mechanics is based on path integration where the quantities sought are 
transition amplitudes between pure states. But the identification of the causal structure of a quantum system requires 

to go beyond such transition amplitudes. In fact, retardation and causality are expressed by means of expectation 
values corresponding some intial state which is not necessarily stationary. Another, related problem of the transition 
amplitude formalism is that Wick theorem holds for time ordered Green-functions only and we have no easy way of 
writing retarded Green-functions in perturbation expansion as a sum of Feynman graphs, an indispensable scheme to 
the intuition. Hence we shall consider expectation values within the GTP scheme, proposed by Schwinger long time 
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ago [6] and has been rediscovered in a number of occasions [7-15]. The main points, used in deriving the effective 
system dynamics can be summarized as follows [16, 17]. An expectation value of an observable O taken at time t can 
be written in the Schrodinger representation as 

(0(t)) =Tr[0e-*(*-*^)'f^/9ie*(*-*^)-f^], (33) 

where pi is the density matrix in the initial state at time ti. Its expectation value and Green- functions can be obtained 
by means of the generator functional 

gtivb"] ^ TrT[e"^ ^'^i '^*[^(*)--'"^(*)°(*)l]pi(r[e"^ ^'^i dt[H(t)+j- (t)o(t)]^^^\ ^^^^ 

given in the Hcisenbcrg representation where j = {j^ is a CTP doublet of sources. W[j] generates the connected 
CTP Green- functions of 0{t). The extension of the single time axis to two, oppositely oriented time axes, corre- 
sponding to the two time evolution operators in Eq. (34) allows us to generalize the Wick theorem for the extended 
time contour and to establish perturbation expansion, based on CTP Feynman graphs. The argument is based on 
the simple observation that Wick theorem follows automatically from a path integral representation of the generator 
functional, 



where the CTP doublet x — (x+,x+) denote a pair of trajectories for each degree of freedom. The name "Closed 
Time Path" refers to the boundary condition x~^{tf) = x~{tf) representing the trace in (34), imposed at an arbitrary 
time tf chosen to be at a later instant than the expectation values taken. The CTP action 

ScTp[x] = S[x+] - tS[x-] + Sbc[x] (36) 

contains the usual action S[x^] of the theory and another contribution, Sbc[x], incorporating the boundary conditions 
in time. Both the initial state of the system, the perturbative vacuum, and the final conditions x~^{tf) = x^{tf) can 
be represented by a quadratic term, S'Bc[a;] = 0{x^). The detailed from of Sbc[x] is not needed as long as the free 
propagators can be constructed in the operator formalism [16] . The CTP generator functional offers a simple method 
to monitor the unitarity of the time evolution. In fact, for unitary time evolution the norm is preserved and W = 0. 
In other words, on has 

W\j,-j] = (37) 

for any Hermitean operator O. 

To find the expectation values we parametrize the sources as = j{l ± k)/2 ± j where k is an arbitrary real 
parameter and write 

m)) = ™ (38) 

which is K-indcpendcnt for unitary time evolution. The source j represents the external environment which generates 
unitary time evolution and guides the system adiabatically to the desired state in which the expectation value is 
sought. The other source, j is a book-keeping device to generate the expectation value and is set to zero after the 
functional derivative in Eq. (38). 

The equations of motion, satisfied by the expectation values can be inferred from the effective actions, obtained by 
means of functional Legendre transformation j -H- A for fixed j, 



r[o] = ^w\j] - J dtj{t)o{t), o{t) = 



SW\j,j] 

The inverse Legendre transform is defined by 



(39) 



?fiWlj] = T[0] + I dxj{t)Oit), - m = (40) 



showing that the expectation value is an extremal of the effective action in the physical case when j = is set. 
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We briefly summarize now the free generator functional and the propagator for a harmonic oscillator and a free, 
real scalar field governed by the actions 



If 1 

Sho[x\ ~ 2 / ^^b^^i^)^ ~ '''^0^^ i^)] — 2^(^o^ ^ 

Sft[4'] = \ j 



%I^oV- (41) 



2 

The generator functional (35) for 0{t) = x{t) and 0{x) = (j){x) can be written as 

■r\J] = j £)[0]et^*o"''^+'-^~^, (42) 
and give 

Wnod] = -^jGoj, 

Wft[J] = -\jdJ. (43) 

The CTP propagator £) of a local observable (j){x) of bosonic exchange statistics can be parameterized by three real 
functions, 

({T[cl>{x)m\) imH^)) \ _ (D^x, y) + iD\x, y) -Df{x, y) + iD\x, y)\ _ ^ 
\ {mm) {T[mm]Y ) ~ \Df{.x, y) + iO^x, y) -D"(x, y) + iD\x, y) J ' 

where D'^{x, y) ~ D"'{y, x) and {x, y) = —D^[y, x) are the near and far field Green-functions and the imaginary part 

is D^{x, y) = D^{y, x). The causal, retarded and advanced propagators are defined as D = + iD', D° = D" ± , 
respectively. Note that the off-diagonal blocks, Wightman- functions, consist of mass-shell amplitudes only. The only 
off-shell contribution of the propagator comes from the near field component and is generated by the time ordering. 
This parameterization can be carried over any two-point function, in particularly self energy, cf. Eq. (57). 
The generator functional (43) assumes the form 

WhoIj] = -jCg - ^jG^j ' ^Gl^j, 

Wft[J] = -JDlJ-'^JD'^J-'-JDiJ. (45) 
when the parameterization (44) is used and the effective actions 

Tho[x\ = ^(x + jG'')G"-i(x + G'^j), 

TftM = ^(<^ + J£>»)D"-i(</, + £)'-J), (46) 



give rise the equations of motion 



-D^'J + kD^'J. (47) 



The propagator is easiest to find in Fourier space, 

G,{t,t') = j ^e-^-(*-*')Go(a;), 



Do{x,y) = I 7K^e-''^<-^-y^Do{k), (48) 
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where 



UJ^ —Lj'f.-\-i€ 



= (-2...(^)e(.0) ^) . (49) 

Wc shall need the inverse of the free propagator what can be obtained by using the regulated expression (5^(w) = 
+ 77^)7r for the Dirac-delta in (49) with the result 

b-^\k) = I?o-i(fc)a + ze(_2Q^^^O) -2e(-fc0)^ ^ ^^^^ 



where the matrix 

/i n \ 

(51) 



1 
-1 

acts on the CTP indices. Note that the transformation e — )■ — e amounts to a complex conjugation of the inverse 
propagator and therefore of the propagator itself, as well. Therefore the transformation e — >• — e preserves G" but 
changes the sign of . 



B. Quantum Harmonic Oscillators 

We use the CTP scheme first for the system of harmonic oscillators whose Lagrangian (1) will be written in 
condensed notation as 

y\ = V +^13 X^G'^Xn - ^ gn^nV (52) 

n n 

without the source term where G~^{uj) = m{u)^ — ui^). The generator functional (35) reads as 
by suppressing Sbc- 

When the dynamics is given in terms of normal modes as in the Lagrangian (13) then the spectral strength (18) 
can be used to find the system propagator, 

= / .02«,„,n) (_,„,pfcji)e,., -""'iti^f '-"') ■ i^) 

The propagator is causal as long as the frequency integral in calculating G{t,t') as in Eqs. (14)-(15) is made before 
the integration over the spectral variable. When the integration is first carried out for the spectral variable then 
sufficiently singular spectral weight, such as given by Eq. (24) give acausal propagation. 

Another way to solve the model which is better suited to deal with interactions is the integration over the environ- 
ment variables to find the influence functional [18], 

= e-^^^y. (55) 

It defines the self energy 

E = alia (56) 



with 



The form 



of the generator functional justifies the use of the dressed propagator 

1 



G 



The propagator, shown in the first line of (49) inserted into Eq. (57) gives 



^jr^^, -27ri5(a;2-u;^)e(-a;)^ 
2 ,a\o,(,,\ 1 



which can be written in terms of the environment spectral strength (16) as 

,2 n2 



fl = /.^n2np(0)(_2^.^P^|^Q^^ 
and straightforward integration of the form (16) gives 

n"(a;) = n"(a;)±n/(w) 

^ ^ mef^D(w2 + f^|,)' 

The inversion of the 2x2 matrix in Eq. (59) is trivial and the result is 



1 



2 



meQoi^D T ?w) 



Go' - m 



C. Quantum field as environment 

Now we turn to the model defined by the action (25), written in a condensed notation as 

S[<t>, y] = lyGo'y + ^ W V + My- 

without the source term where the free system inverse propagator is given by Eq. (9), the free 
propagator is 

2„2 



and 



A{k,Lj) = —g6 ( w I . 



The generator functional (35) assumes the form 
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TABLE I: The renormalized parameters of the retarded Green-function 



_d rwd/m c^d/i^o 7d_ 

11 1 1 



3 1 1-^^ 



11 1 



1 f ('i£n£r s V ■'(-> I 1 



TV vq 47r 

mi 



where 

Go He.) = {us' - + ie (^_^^^^^ -2e{-us)\^ ^^^^ 

denotes the free system propagator and the inverse environment propagator Dq^ is given by Eq. (50). The effective 
theory for y can easily be obtained by integrating out ^, 

This result leads to Eqs. (59)- (56) with ft = A*'' DA. 
The actual form of the self energy 



n(u;) = g^c^ / -—: ^ -,7 2 ^%]\^r^f ^ 1 • (70) 



yields the retarded, advanced and imaginary parts, 



n»(w) = 52c2[5i/^(^)±isign(a;)9/d(w)], 

W{u) = g^c^QIdicj), (71) 



with 

d'^k 1 

'\k\ 

A being an UV cutoff. The loop integral for d = 2, 3, 4 and 5 turns out to be 

i 



f df'k 1 



h 



2c\uj\ ' 
1 



47rc2^''A2c2 4c2' 
A i\uj\ 

A3 Ac^2 



727r3c2 247r3c4 247r2c5 ' 



(73) 



in the case of a gapless environment, = 0. Note that the infinitesimal imaginary part of the free inverse propagator 
(68) can be omitted in the denominator of the dressed propagator, (59) because the self energy represents a similar, 

finite contribution. 

Insertion of Eqs. (71) and (73) into Eq. (63) represents our final result for the propagator in odd dimensions, 

1 
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and in particular , 

G" = 



...2/ 9 ,2^,2 I ^2^2,,,4/u^X2d-4- ^'°> 



The length scale £ was introduced after Eq. (25) to parameterize the coupling constant g and to^, and 7^ are listed 
in Table 1 where the UV cutoff is given in terms of the minimal distance Tq = 27r/A. The even dimension, d = 2 we 
find logarithmic terms, 



1 



2 / „ o III? , ,.,2 \ tsign(ai)g^ ' 

^ A 



m (a;2 _ ^^2 _ g 



^2 = : -2 7. (76) 

_1_ 9l 



Aw? {u 



47r (?A^ ; ' 4 



The parameter £ can be interpreted as the characteristic length of interactions. In fact, the observation of the system 
with a period r in time leads to a self energy proportional to {£/ctY~'^ in odd dimensions. The system-environment 
interaction is UV or IR dominated for d > 2 or d < 2, respectively. The renormalized parameters of Table I depend 
on I/tq, the interaction length counted in units of the minimal length and 1/ctq where Tq denotes the time scale of 
the sampling of the environment by the system. The UV renormalization is weak in the super-Ohmic regime, d < 3 
but becomes strong and power like for d > 3. 

Note the equivalence of the classical and quantum effective theory for quadratic systems, the identity of the retarded 
self energies, Eq. (30) and If" + 11^, given by Eq. (71), respectively. 

The Green- functions (74)- (76) are acausal except in the Ohmic case, d = 2> where the only compromise between 
causality and dissipation, an equation of motion which is quadratic in the frequency with linear damping term is 
realized. The normal frequencies are 

w± = (± Vl - a? - ia)ujs (77) 



where a = ^/87ry 1 — 7;;^. The theory with ^ < Stt displays a quantum phase transition at a = 1, the system 

becomes overdampcd for a > 1. In other odd dimensions there are poles on the wrong side of the real frequency axis 
and causality is lost. One expects dissipation and imaginary, odd terms in the frequency in the denominator in any 
odd spatial dimensions apart of three. It remains to be seen what is left from the causal structure in even dimensions. 

When the environment has a gap in its excitation spectrum, H ^ Q then the loop integral (72) and no dissipation 
or acausality is observed at low enough frequencies, [a;| < fic^ /h. Dissipation is absent in a similar manner in case of 
a finite system because the energy may flow back to the system and establish an equilibrium. 



IV. NON-PERTURBATIVE RADIATION DAMPING 



We have derived the Green-functions for the system coordinate y by considering classical and quantum effective 
theory obtained by eliminating the environment. The system coordinate y = G^j satisfies the equation of motion 

j = {GT'y. (78) 

What initial conditions do we have to provide to make the solution unique? The conditions given by Eqs. (3) may 
not be sufficient because the effective equation of motion, Eq. (78), usually contains higher order derivatives in time. 
To find the missing conditions let us return to the case where the environment is made up by a finite, N number 
of oscillator and the system coordinate is given by Eq.(7). Initial conditions for the environment variables appear 
through Xhnit') in the solution. Therefore the initial conditions (3) together with the other 2N conditions for the 
environment can be adjusted to fix the system coordinates and its first 2A^-|- 1 derivatives at the initial time. We have 
the opposite problem when the effective equation of motion is known only without the dynamics of the environment. 
In this case the informations in the initial conditions required beyond (3) actually identify the initial conditions for 
the environment. 
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Such a transmutation of the initial conditions is more involved in an infinite system. In our field theory models 
defined in an infinite volume a part of the environment excitations generated by the system remain localized in space 
and another part moves away and propagates to infinity. Such a decomposition of the induced field is given by 
the near and the far field Green- functions. In fact, the ofi^-diagonal CTP blocks of the propagator are on-shell as 
mentioned after Eq. (44) and describe the radiation field. The off-shell amplitude of the near field Green-function 
produces localized response. Initial conditions are introduced for the system and the radiation field only. Therefore 
the additional initial conditions needed to fix the solution of the effective equation of motion (78) are belonging to 
the radiation field, described by , shown in the second equation of Eqs. (75). Hence the higher order terms in 
the effective equations of motion represent the interaction of the system with its own radiation field. The near field 
renormalizes only the parameters of the free equation of motion. 

The impact of the induced field on the motion of a point charge is well known in classical electrodynamics, it is 
the Abraham-Lorentz force [19]. Its actual form is analogous to the scalar field in d = 5 because the vector-current 
minimal coupling contains a space-time derivative. We need a single additional initial condition in this case, it is is 
enough to parameterize the infiuence of the radiation field on the initial condition of a point charge. The absence of 
higher than third order derivatives with respect of the time in the Abraham-Lorentz force can be understood easily by 
following an expansion in the retardation [20] where the higher order contributions are suppressed in the point charge 
limit. The role of the initial conditions mentioned above suggests that the run-away solutions should be avoided in 
electrodynamics by the use of appropriate intial conditions for the particle trajectory. 

Perturbation expansion is usually employed in the interactions rather than in the retardation. Note that this is not 
a useful scheme for higher order derivative terms in the linearized equation of motion. When these terms arc treated 
as perturbation then the unperturbed equation of motion remains second order and no additional initial conditions 
are needed to identify the solution but we have to specify completely the initial radiation field. These latter explains 
the absence of perturbative run-away solutions. The dissipation is a non-perturbative phenomenon, as well, because 
the perturbation expansion breaks down after long enough time and it can not give account of irreversibility, the 
impossibility of recovering the radiated energy loss of the charge from an infinite volume. Acausality is another 
non-perturbative effect because the location of the poles of the retarded Green-function on the complex energy plan 
is not an analytical function of the coefficients of the higher order derivative terms. 



V. COARSE GRAINING IN TIME 



The examples above show that there are two different realizations of the dynamics, depending on the order of 
integrations over the frequency and the normal modes in Green- functions. We now locate the source of this ambiguity 
in mathematical and physical terms. 

It is easy to find the mathematical source of the ambiguity. Fubini's theorem assures that the order of integration 
is arbitrary as long as the integrand is a continuous function of the integral variables. But the Fourier transform of a 
Green-function of continuous spectrum is singular in a region called mass-shell and the result, obtained by performing 
the integrations in different orders may differ. There are two infrared cutoffs in the dynamics, one is the spatial size 
L which is related to the number of particles N of fixed density and the time of observation T. The removal of these 
cutoffs, the limits 



and 



lim lim (O) (79) 

N->-oo T->-oo 



lim lim (O) (80) 

T-foo JV->-oo 



may differ. But what is then the correct order? 

Mathematical ambiguities always point to a choice in the preparation or in the observation of the system. It is the 
shape of the space-time region where the solution is sought which introduces this freedom in the present case. One has 
to integrate first over the energy-momentum components with continuous spectrum and sum of the discrete spectrum 
next. When the size of the space-time region diverges in several directions then the component corresponding to 
the faster growing size or to the faster decreasing level spacing of the spectrum is integrated before the others. In 
other words, the order of integration over the continuous energy-momentum spectrum of a theory defined in infinite 
space-time depends on the relative speeds of the limits T — >• oo and L — >■ oo. 

One is usually interested in the dynamics of a system enclosed in a large but finite quantization box. According to 
the standard rule of quantum field theory one carries out the frequency integrals before the momentum integration. 
But observations are made in finite length of time with an unavoidable uncertainty in measuring frequency or energy. 
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There are two circumstances when this hmitation is important. One is when we encounter unusually slow degrees 
of freedom. An example of this case is spontaneous symmetry breaking where the time dependence of the order 
parameter experiences critical slowing down in the vicinity of a second order phase transition and any finite restriction 
on the observation time leads to the recording of apparently constant order parameter with symmetry breaking value. 
Irreversibility is actually related to a dynamical symmetry breakdown and appears as a phase transition. Initial 
conditions imposed on a dynamical system do not lead to the loss of time reversal invariance because they leave the 
equations of motion unaffected. But the elimination of some dynamical degrees of freedom, the environment, generates 
nonlocal effects in time. Hence the initial conditions of the environment may break the symmetry under time reversal 
in the effective system dynamics. The complication in establishing such a spontaneous symmetry breaking arises from 
the difficulty of finding a local order parameter, the entropy being a non-local quantity. 

Another, more phenomenological consideration which makes the imprecise measurement an important issue is 
related to incomplete information. A loss of information about the dynamics can be viewed as a coarse graining either 
in space or in time. It is well known that a coarse graining in space, the ignorance of spatial correlations or degrees 
of freedom may lead to irreversibility when the system becomes infinite and the slowing down of a certain part of the 
dynamics manifests itself as a diverging Poincare recursion time. Coarse graining can be imagined in time, as well. 
Informations about the environment are encoded in the complexity of the effective system dynamics and we need 
precise observations to regain them. One has to resolve all environment oscillator frequency to restore cornplcitcly 
the environment dynamics. This can be achieved by observations of finite length as long as the spectrum is discrete. 
But such observations leave infinitely many environment degrees of freedom unresolved when the spectrum has an 
accumulation point. Therefore a condensation point of the spectrum at the ground state may imply irreversibility. 

To understand this point better we implement a long time cutoff in the model defined by the Lagrangian (1) by 
assuming that one observes the product 



of the system coordinate and an IR cutoff function satisfying the conditions c{t) ~ 1 for t <C T and c{t) ~ when 
t^T. Let us suppose that the external source is localized in time, j{t) = 6{t). Then the observation of yobs{t) allows 
us to construct the Fourier transform 



We shall attempt to draw conclusion about time reversal invariance and causality by extracting the environment 
spectrum from (82) exclusively. 

When N is finite and we have unlimited time of observation at our disposal, T = oo then c{lu) = 1 and we can 
resolve the discrete spectrum. As the time of observation starts to be limited the observed spectrum spreads. In fact, 
let us use the Green- function (19), 



y{t) Vobs 



it) = c{t)y{t) 



(81) 




(82) 




(83) 



and shift the integral variables, oj' ^ co' + co, fl ^ fl + aco' 



''. The result is the expression 




(84) 



including the apparent spread spectral functions 




(85) 



Let us choose for the sake of definiteness the IR time cutoff 




(86) 



The spectral functions 




(87) 
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(a) (b) (c) 

FIG. 1: The apparent spectral function pj^' for m = 1, uij = ujo/j, N = 20, (a):T = 2000, (b):r = 700, (c):r = 100 as functions 

of Uj/uiQ. 

are plotted the simple spectrum Uj — ojo/j, 3 — 2, . . . and mixing Aqj = 1/N for different N in Fig. 1. For long 
enough observation time, T — >■ cx), the discrete peaks can be resolved in Figure 1(a). But whatever long observation 
time is allowed, a condensation point in the spectrum always yields unresolved peaks. In fact, the peaks at frequencies 
ij and u}' can be resolved when r^l/|a) — a)'|, a condition not satisfied by any pair of peaks by finite T in case of 
a condensation point. This phenomenon is already visible as the slight increase of the local minima of the spectral 
function in Fig. 1(a) around w/wq ~ 1.05. The overlap of the peaks is more visible for shorter observation time as 
shown in Fig. 1(b). Finally, for an even shorter observation time, shown in Fig. 1(c), there is hardly any evidence 
left from the discrete nature of the spectrum. 

The choice between the schemes (79) and (80) is now simple. If there is time to resolve all peak, for instance 
for discrete spectrum, then one should apply the prescription (79). If the observation time is too short to recognize 
every peak then the procedure (80) is justified. But no finite amount of time is enough to resolve all peak when the 
spectrum has a condensation point. More precisely, finite time observations resolve finite number of degrees of freedom 
and the rest is dumped into the continuous part of the spectrum. Such a loss of information leads to an apparent 
breakdown of time reversal invariance when infinitely many unresolved peaks suck the energy from the resolved part 
of the dynamics without having time to return it. 

What is the time needed to discover the spontaneously generated time arrow? The shift of poles due to the spread 
of the discrete spectrum generated by a large IR cutoff T is 0{T~^). Irreversibility appears through the finite life 
time of excitations which is 0{T) in this case, betraying its artificial origin. When there is a part of the spectrum 
which remains unresolved then its analytic continuation may generate poles at distance larger than 0(T~^) from the 
origin and may shift poles by more than 0{T~^). Consequently, the life time of elementary excitations can be shorter 
than 0(T) and irreversibility may become a genuine, cutoff independent phenomenon. 

It is instructive to compare explicit and spontaneous breakdown of the time reversal invariance within our Gaussian 
model. An explicit irreversibility, achieved by adding a finite, non-Hermitean term to the Hamiltonian leads to 
complex action. In the case of spontaneous symmetry breaking of time reversal invariance the action is real. Hence 
the equation of motion is a real function of iu) and poles of the Green functions with non-vanishing real part appear 
in pairs, uj± = iwi + iuj2- The term with alternating sign, ±uJi is the remnant of the formal time reversal invariance 
of the original, elementary dynamics and the imaginary part, iu}2 refiects irreversibility. 

Explicit breakdown of time reversal invariance leads to non-unitary time dependence in quantum models. But the 
spontaneous breakdown of the time reversal invariance by soft modes preserves formal unitary time evolution as long 
as all degrees of freedom are followed. In fact, the form (45) of the generator functional is always valid and gives 
W = for arbitrary j when j — is imposed. The system state "leaks" to the environment, indicated by the complex 
poles of the system Green-functions but unitarity is maintained on the level of the whole system. Naturally one needs 
infinitely long time to confirm unitarity by observations, the detailed time evolution appears non-unitary in any finite 
time interval. 

Causality can be tested by monitoring the response of a source localized in time before its action. The translation 
invariance in time is broken by the initial conditions and such a test has to be performed by acting with the source at 
different times, rendering causality a non-local feature in frequency space. A more formal way of seeing this is to note 
that the location of the poles of the free Green-functions is a non-analytic function of the coefficients of the higher 
order time derivative terms in the free Lagrangian. The uniqueness of the analytic continuation of the propagator on 
the complex energy plan protects causality when the spectrum is discrete. But condensation points or a continuous 
part in the spectrum may shift poles to the nonphysical sheet. In more physical terms, the abundant soft modes may 
render the causality test unreliable. In fact, one attempts in this test to determine the time of an external source by 
observing the system coordinate and a restriction on the observation time may make the reconstruction of the time 
of the external perturbation difficult when the dominant modes are almost time- independent. 
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VI. SUMMARY 

The origin of irreversibility and acausality is identified in this work in the framework of simple harmonic systems. 
The detailed equations of motion of the model are causal and appciar formally invariant under time reversal but 
the effective dynamics of a single degree of freedom may be irreversible and acausal. In other words, time reversal 
invariance and causality might be lost when infinitely many reversible and causal degrees of freedom are properly 
cooperating. Such a loss of symmetry which is present in the elementary equations of motion is a reminiscent of 
spontaneous symmetry breaking. The initial conditions imposed on the environment generate a dynamical breakdown 
of time reversal invariance which may persist in time in the when the environment modes are sufficiently soft. 

Both irreversibility and acausality can be related to the insufficient frequency resolution in the infrared. This issue is 
negligible for systems with discrete frequency spectrum but becomes important when the spectrum has a condensation 
point at the ground state. The limited freqiiency resohition leads to extrapolations which suggest the spread of the 
spectrum lines and the presence of dissipative forces. Verification of causality, the temporal order of the action of an 
external source and the appearance of its consequences is rendered unreliable by the unresolved, continuous part of 
the spectrum consisting of almost time-independcait modes. 

The breakdown of time reversal invariance, the appearance of complex poles in the Green-functions is usually 
followed by acausal behavior in lacking of a specific mechanism to protect the unphysical sheet against the intrusion 
of poles. Within the class of the models considered the only exception was the massless scalar field environment in 
three spatial dimensions where irreversibility is reflected by the lowest odd power of the frequency in the equation of 
motion and the pole is happened to stay on the physical sheet. 

We do not attempt to clarify the origin of the arrow of time in this work. But it can be seen in our calculation 
that the time arrow for the system is set by the initial conditions imposed on the environment. This is a mechanism 
which, in a manner similar to the sensitivity of an unstable equilibrium position against fluctuations may spread a 
time arrow, set for some degrees of freedom by a sufficiently restrictive initial condition over the rest of the system 
irrespectively whatever weakly may they be coupled. 
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